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We present an efficient scheme for the controlled generation of pure two-qubit states possessing 
any desired degree of entanglement and a prescribed symmetry in two cavity QED based systems, 
namely, cold trapped ions and flying atoms. This is achieved via on-resonance ion/atom-cavity 
couplings which are time-dependent and asymmetric, leading to a trapping vacuum state condition 
which does not arise for identical couplings. A duality in the role of the coupling ratio yields states 
with a given concurrence but opposing symmetries. The experimental feasibility of the proposed 
scheme is also discussed. 



Entanglement is one of the most subtle but striking 
quantum phenomena. In the quantum information pro- 
cessing (QIP) field, it is considered as a physical resource 
which can be quantified and manipulated. The power 
of this quantum resource has been demonstrated in var- 
ious applications such as state teleportation, quantum 
cryptography and quantum dense coding A crucial 
ingredient for any QIP protocol is the controlled and ac- 
curate generation of entangled multi-qubit states. A fast 
and deterministic generation process is highly desirable 
in order to minimize decoherence effects. 

Strongly coupled ion/atom- photon systems are prime 
candidates for QIP Many schemes have been pro- 
posed for entangling atomic qubits using cavity QED, 
with some even being implemented experimentally 2, 
0, 0]. Zheng et al. proposed entanglement via a vir- 
tually excited cavity mode Q which was later experi- 
mentally realized |4j. Specifically, two atoms crossing a 
non-resonant cavity are entangled by coherent energy ex- 
change. However, resonant cavity QED is likely to offer 
faster entanglement schemes. Plenio et al. considered 
singlet-state generation through continuous observation 
of the cavity field in resonance with two two-level atoms 
0. The resulting entanglement generation is, however, 
comparatively slow and probabilistic, requiring efficient 
active measurements. Most theoretical studies on cavity- 
assisted entanglement assume identical qubit-cavity cou- 
plings. Nevertheless, it is known that spatial and/or tem- 
poral variations of the ion/atom- field coupling arise ex- 
perimentally HIE 0- Zhou et al. recently considered 
entanglement within a non-symmetric two-atom system 
coupled to a thermal cavity state |8| , and found that the 
entanglement is maximum for a particular static asym- 
metry in the coupling. However the full consequences of 
having dynamical qubit-cavity couplings, are unknown. 
Indeed, most existing proposals are aimed at the gener- 
ation of maximally entangled two-qubit states without 
explicit reference to the resulting state's symmetry. This 
motivates us to explore systems in which the interplay 
between entanglement and state symmetry can be stud- 



ied directly, with a view to extracting these states for 
QIP and to gain insight into physical mechanisms for en- 
tangling large numbers of quantum subsystems. 

Here we propose a scheme to generate and control two- 
qubit entangled states possessing a definite symmetry 
which can be chosen a priori. Our results are obtained 
using a generalized Dicke model with time- dependent 
qubit-cavity couplings, thereby mimicking actual experi- 
mental conditions 01(3,13 • The resulting diversity of cou- 
pling profiles and timescales, offers a natural mechanism 
for manipulating the collective qubit dynamics and hence 
controlling the entanglement. We show that it is the geo- 
metric mean of the individual couplings which dominates 
the entanglement's dynamical evolution. Exploiting a 
particular trapping vacuum condition in which the 
cavity-qubit state becomes separable but the two-qubit 
state remains entangled, we show how full control of the 
quantum correlations can be achieved. We then demon- 
strate the deterministic generation of a pure two-qubit 
state with any degree of entanglement and a prescribed 
symmetry. Our scheme has the following advantages: (i) 
since it is resonance-based, the entanglement generation 
requires a significatively shorter operation time than both 
the ion/atomic dipole decay time and the radiation de- 
cay time; (ii) the resulting generation of pure states with 
an arbitrary degree of entanglement, is deterministic and 
does not require final state projection - hence it is accu- 
rate; (iii) since the system is non-symmetric with respect 
to qubit exchange, the final-state symmetry can be con- 
trolled by choosing which qubit should initially be in its 
excited state. This final-state symmetry can therefore be 
used to encode any quantum information present in the 
initial atomic state. 

We investigate the dynamical evolution of the concur- 
rence C(p n . „.) which measures the degree of two-qubit en- 
tanglement |lOj, and the two-qubit correlation function 
{ <7 2 <J i) which measures the degree of qubit exchange 
symmetry. For a singlet state (c^~err) = —1/2, while 
for a maximally-entangled triplet state (c^ct^) = 1/2. 
Hence negative values for (c^cr-j - ) are associated with 
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an anti-symmetric-like behavior while positive values are 
related to a symmetric-like behavior. Furthermore, we 
also investigate the qubit-field correlations {afa) and 
find relations between these quantum correlations and 
the qubit-field entanglement. 

We consider two possible scenarios: (i) cold trapped 
ions inside a high-finesse cavity where the Lamb-Dickc 
localization limit is assumed and (ii) flying two-level 
atoms crossing a cavity ^2|. I n both cases, given that the 
atomic transitions are on resonance with the single-mode 
cavity field, the Hamiltonian in the interaction picture 
and rotating- wave approximation becomes {% = 1) 

Hi = fi(U,t, r t ){aVr + a+a} (1) 

2=1,2 

where at = |e<> <£T< 1 , ^ = |ffi)( e *l with |e<) and \g^ 
(i = 1,2) being the excited and ground states of the 
i'th qubit. Here a) and a are respectively the creation 
and annihilation operators for the cavity photons. The 
time-dependent coupling of the cavity field with the i'th 
qubit, which is injected at U and interacts during a time 
Tj, is given by a time- window function, 

fi(U,t,Ti) = [Q(t - n ) - Q(t -n- U)Mt) (2) 

where ji(t) is the time-dependent qubit-field coupling 
strength. We focus here on the situation in which the 
qubits interact simultaneously with the cavity mode such 
that ti — t 2 — and t\ = t 2 = t. 

In the case of identical time-dependent couplings 
fi(ti,t,Ti) = f 2 (t 2 ,t,T 2 ) — f(0,t,r), the Hamiltonian 
is symmetric with respect to qubit exchange and hence 
commutes with the pseudo-spin operator J 2 = J 2 + 
Jy + J 2 where J x ,y,z are the usual angular momen- 
tum operators. The dynamical evolution of the state 
does not mix different J— sectors. Hj commutes with 
itself at different times as well as with the operator 
V = 2'i a ^ cr r + cr i +a }- The combined qubit-field sys- 

tem follows a unitary time evolution generated by the op- 
erator U = exp(—i(p(t)V), where 4>(t) = f Q j(t')dt' repre- 
sents the net effect of the time-dependence in the qubit- 
field interaction. For different but time-independent cou- 
plings, i.e. 71 ^ 72, the Hamiltonian has two impor- 
tant properties. First, it no longer commutes with J 2 
and second, it is non-symmetric with respect to qubit- 
exchange. We will show how these features can be ex- 
ploited to control the total pseudo-spin and hence the 
symmetry of the two-qubit state. The time-evolution 
of the system is more complex for the case of different 
but time-dependent couplings, since generally the Hamil- 
tonian does not commute with itself at different times. 
Hence a history-dependent dynamics arises, governed by 
U = Texp(—i f Q Hj(t')dt r ) where T denotes the time- 
ordering operator. In all these cases, the number of ex- 
citations N = a) a + X)j=i 2 °i IS a conserved quan- 
tity. This implies separable dynamics within subspaces 
having a prescribed eigenvalue N of M . For N = 1 a 
basis is given by {\e 1 ,g 2 , 0), \gt, e 2 , 0), \gt, g 2 , 1)}, while 



for N > 2 a basis is given by {|ei, e 2 , N — 2), \ei,g 2 , N — 
1), \gi,e 2 , N—l), \gi,g 2 , N}}. The third label denotes the 
number of photons. For simplicity, we will restrict our- 
selves to consider |\&(0)) = |ei,<72>0) as the initial state so 
that the system's dynamical evolution is confined to the 
single excitation subspace. The whole system's quantum 
state at any time can be expressed as 

|*(t)> = oi (t) |ei, fl2 , 0) + 02 (t) \gi , e 2 , 0) + a 3 (t) \gi , g 2 , 1) (3) 

Clearly the cavity mode acts as a third qubit, 
hence the single-qubit-field entanglement can be quan- 
tified by the concurrence qubiti^-field as well. We 
have found simple relations between these quanti- 
ties: C(p a , a )(t) - 2\(a+(t)ai(t))\ = 2\at(t)a 2 (t)\, 
C{ P ai.f)(t) = 2\(a+{t)a(t))\ = 2K(i)a 3 (i)| and 
C( P a2j)(t) - 2\(a+{t)a(t))\ = 2\a* 2 (t)a 3 (t)\. We con- 
sider the situation in which 71 (t) = r"/ 2 (t), where r is a 
constant. The unitary evolution is given by 

oti(t) = l + ra 2 (t) 

a 2 (t) = -2aSin 2 [6»(<)/2] (4) 
a 3 (t) = -iy/faSm[0(t)] 

where 9{t) = J* uj(t')dt' is the effective vacuum Rabi an- 
gle. The time-dependent frequency of the collective qubit 
mode coupled to the cavity field is given by w 2 (i) = 
7i W 2 + 72W 2 , and a = 7 1 (i) 72 (i)/c i ; 2 (i) = r/(l + r 2 ) 
denotes the relative geometric mean of the couplings. 

To help understand the physics for non-symmetric cou- 
plings, we express ^(t)) as 

\V(t))=A(t)\$-,0) + B(t)\$ + ,0)+a 3 (t)\ 9l ,g 2 ,l) (5) 

where |$±) - [| ei , 52 ) ± \gi,e 2 )]/V2, A{t) = [a 2 {t) - 
ai(t)]/V2, and B{t) = [a 2 (t) + ai(f)]/v^. From Eqs. 
and JSJ it is clear that when 0(r*) = n the qubit-field 
state is separable and is given by |\£(t*)) = [A(r*)|(f> _ ) + 
S(r*)|$+}] ® |0). Therefore, at this special time a trap- 
ping vacuum state condition holds, i.e. the cavity photon 
number is unchanged. For identical couplings, i.e. r = 1, 
the state |$~, 0) is an eigenstate of Hj(t), thus the coef- 
ficient A(t) remains constant in time, i.e. A(t) = A(0) = 
-1/V2. Additionally for t = r*, B(t*) = l/y/2 such that 
the system's state becomes fully separable, i.e. |$(r*)) = 
— |gi,e2,0), hence there is no entanglement in the qubit 
subsystem nor between the qubits and the field. Figure 
^ shows the dynamical evolution of the two-qubit con- 
currence and the qubit-field concurrences, for the case of 
identical time-dependent couplings. The initially excited 
qubit becomes entangled with the cavity field faster than 
the initially non-excited qubit. Nevertheless, there is a 
symmetry in the dynamics of individual qubit-field en- 
tanglements - both qubits become identically entangled 
with the field but at different times. At t = r*/2 when 
the individual qubit-field entanglements are identical, i.e. 
C(paij) — C(pa2,f), the two-qubit concurrence reaches 
its maximum value C(p Qia )(r*/2) = 1/2 corresponding to 
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the state |*(r*/2)) = [-\$-,Q)-i\g 1 ,g 2 ,l)]/V2- Hence 
for r = 1, two-qubit entanglement requires the existence 
of qubit-field entanglement as well. 

With different couplings, the symmetry with respect 
to qubit exchange is broken. Hence the coefficients 
A(t) and B(t) can be fully controlled by choosing ade- 
quate values for the coupling ratio r. In particular at 
t = t*, any allowed superposition of |$~) and |$ + ) is 
obtainable in such a way that any prescribed two-qubit 
entangled state can be achieved. Figure shows the 
two-qubit concurrence as a function of r. Two impor- 
tant physical consequences follow from Figure [3 First, 
the symmetry of the two-qubit state is controlled by 
the parameter r. If the initially excited qubit is in- 
teracting with the cavity through the weakest coupling 
(r < 1), the two-qubit state is in the anti-symmetric sec- 
tor ((<rt( T *) (7 i( T *)) < (FigEJa)); otherwise it is in the 
symmetric sector ((cr+ (r*)oT (r*)) > (FigEJb)). Sec- 
ond, when the trapping vacuum condition is met (t = r*), 
it is always possible to get two different matter states 
with identical concurrences: one has a positive corre- 
lation (r > 1) while the other has a negative correla- 
tion (1/r). However, these two two-qubit states have the 
same value of the relative geometric mean of the cou- 
plings a = r/(l + r 2 ). Note that a < 1/2, where the 
equality holds for the identical coupling case. Therefore 
for any given value of a, two-qubit states with identical 
concurrence but different symmetries can be found. In 
particular, r — y/2 + 1 and r = V2 — 1 yield the same 
a = for which C(p a . a ) = 1. Hence two maximally 
entangled states can be generated from a single excita- 
tion: the symmetric one for the former value of r and 
the fully anti-symmetric one (singlet state) for the latter 
value of r. Thus, a is the relevant control parameter for 
generating a two-qubit state with any prescribed degree 
of entanglement. In FigE^c) the two-qubit concurrence is 
shown as a function of a. It is a non-monotonic function 
of a as given by C{p a , a ){ct) = 4ay/l — Act 2 . 

Figures and 01 show the dynamics of formation of 
maximally- entangled qubit states with different time- 
dependent couplings. The two-qubit concurrence and 
qubit-field concurrences are shown as a function of time 
for r = V2 - 1 (FigEJ and r = y/2 + 1 (PiggJ). In 
the formation of the singlet, the state remains within 
a negative region of symmetry (Fig|3Jb)), while for the 
r > 1 case the degree of symmetry evolves from negative 
to positive (Fig0Jb)). In both cases, the maximal two- 
qubit concurrence is reached at the moment when both 
individual qubit-field concurrences vanish. It is worth 
noting that the time formation for the anti-symmetric 
state is always larger than the time required for ob- 
taining the symmetric one. For example, for the time- 
dependent couplings depicted in FigOJa) and given by 
7i(t) = ImaxSin 2 (irt / V*) , the formation time for the anti- 
symmetric state is t* s = (2n / '^™ lax } ^fajr and the time 
required to obtain the corresponding symmetric entan- 
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FIG. 1: Two-qubit concurrence C(p a , a ), qubiti-field concur- 
rence C(p a i,f) and qubit2-field concurrence C(p a 2,f) as a 
function of time, for identical time-dependent couplings r = 1. 
The initial state is |ei,g2,0). Inset (a): Coupling strength as 
a function of time. Inset (b): Time variation of the two-qubit 
correlation (cjcri~). 
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FIG. 2: Two-qubit concurrence C(p a ,a) as a function of the 
coupling ratio r, for t — r* and initial state |ei,(?2,Q). Inset 
(a): Two-qubit correlation function (cj (t*)oT(t*)) f° r r < 1 
and inset (b): (aj(r*)<r^(r*)} for r > 1. Inset (c): C(p a ,a) 
as a function of the relative geometric mean of the couplings 
a. 

gled state is r* = rr* s (N.B. r < 1). Although a similar 
dynamics can also be observed with time-independent 
couplings by keeping the same coupling ratio r (trapping 
condition time t* = 7t/cj), it is interesting to note that 
time-varying couplings yield to a switching behavior of 
the two-qubit concurrence dynamics which could be im- 
portant for minimizing decoherence effects (Fig|3{c)). 

In the experimental set-up described in Ref. p], scan- 
ning of an optical field is performed by translating the 
cavity assembly relative to the trapped ions' position. 
Therefore, deterministic control of the spatial/temporal 



4 




FIG. 3: Dynamics of formation of the singlet state |$~,0). 
Two-qubit concurrence C(p a ,a), qubiti -field concurrence 
C(paij) and qubit2-field concurrence C(p a 2j) as a function 
of time, for r = y/2 — 1 . The initial state is |ei,g2,0). Inset 
(a): Coupling strengths as a function of time. Inset (b): Time 
variation of the two-qubit correlation function (crjcrj"). In- 
set (c): C(p a ,a) for different but time-independent couplings 
(dashed-dotted line) and for the time-dependent couplings as 
depicted in inset (a) (solid line). In both cases r = >/2 - 1 
and the time r* is identical. 




FIG. 4: Dynamics of formation of a triplet maximally en- 
tangled state |$ + ,0). G[p a ,a), C{p a \,f), C(p a 2,f), and initial 
state as in Fig. 3, but now r = v2 + 1. Inset (a): Coupling 
strengths as a function of time. Inset (b): Time variation of 
the two-qubit correlation {a^ co- 



variation of the ion-field coupling is feasible. This ex- 
perimental situation is likely to generate couplings of 
order j m ax — 2tt x 3.9MHz for calcium ions (D 3 / 2 to 
P1/2) with a dipole decay rate r = 2tt x 1.7MHz and a 
cavity-photon decay rate n = 2ir x 0.24MHz 13] . The 
interaction time required to generate a symmetric maxi- 
mally entangled state (r = y/2 + 1) is r* = 10~ 7 s, which 



is one order of magnitude shorter than the photon de- 
cay time T r — 1.3 x 10~ 6 s. We note that the required 
time to obtain an entangled state with C(p a , a ) = 0.48 
(a = 0.12), in the symmetric region (r = 8), is even 
shorter: r* = 3 x 10~ 8 s. Our proposal could also be 
realized with flying atoms sent simultaneously through a 
resonant cavity, in such a way that they follow different 
paths hence yielding different temporal coupling profiles. 
Off-resonance coherent control of the collision of two Ry- 
dberg atoms has been experimentally implemented |4J. 
For these atomic Rydberg states, with principal quan- 
tum numbers n = 51 (|e)) and n — 50 (| .<?)), interacting 
with a cavity of quality factor Q — 7 x 10 7 and corre- 
sponding photon lifetime T r = 2 x 10~ 4 s, couplings of 
order 7 maa: = 2ir x 50KHz were achieved. In this experi- 
mental scenario we estimate t* = 10 _5 s for r — \/2 + 1 
and t* = 2 x 10 _6 s for r = 8. These times are obvi- 
ously much shorter than the photon decay time. The 
flying-atoms scheme opens up the intriguing possibility 
of producing a stream of atomic pairs, where each one 
has a different entanglement and/or symmetry. 

In summary, we have proposed a scheme to generate 
pure two-qubit states with an arbitrarily prescribed de- 
gree of entanglement and symmetry. Unlike previous 
cavity-based schemes, the present one takes advantage 
of spatial and temporal variations in the ion/atom-field 
coupling. It is achievable within very short operation 
times since the qubits are on resonance with the cavity 
field. In particular, our scheme is realizable with present 
experimental methods and hence opens up the prospect 
of real-time engineering of multi-qubit entanglement in 
asymmetric time- varying ion/atom-cavity systems. 
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and ORS (AOC), DTI-LINK (NFJ), and COLCIENCIAS 
project 1204-05-13614 (LQ). 
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